Colour Screening, Quark Propagation 
in Nuclear Matter and 
the Broadening of the Momentum Distribution 

of Drell-Yan Pairs 

J. Dolejsf^'^, J. Hiifner^, B.Z.Kopeliovich^'^ 

^ Inst, of Theoretical Physics, University of Heidelberg 
Philosophenweg 19, D-6900 Heidelberg 

^Dept. of Nuclear Physics 
Fac. of Mathematics and Physics, Charles University 
V Holesovickach 2, CS-18000 Praha 8 

^ Laboratory of Nuclear Problems 
Joint Institute for Nuclear Research 
Head Post Office, P.O. Box 79 
101000 Moscow, Russia 

Abstract 

We calculate the broadening of the transverse momentum distribution of 
a quark propagating through nuclear matter. Colour screening plays a fun- 
damental role in that it cuts off quark- nucleon interactions with soft gluons. 
The mean transverse momentum of the quark acquired along its trajectory, 
observed via Drell-Yan pairs, is related to it the ratio of the total inelastic 
meson-nucleon cross section it to the meson mean squared radius. Parameter- 
free calculations agree with the data. 



1 



Nuclei may be used as micro-laboratories for the observation of processes which 
cannot be studied otherwise, e.g. strong interaction dynamics at short-time intervals 
before or after an interaction. This information is completely lost in the scattering 
on a nucleon. Nuclei provide a unique opportunity to study properties of hadronic 
systems at early stages of their formation. Propagation of a quark (antiquark) 
through nuclear matter is one of the examples. The momentum distribution of such 
a quark after the multiple interaction with the target nucleons can be measured in 
a Drell-Yan process of lepton-antilepton pair production. Since the lepton pair does 
not interact on its way out of the nucleus, it carries the undistorted information 
about the interactions of the quark. 

In this paper, we study the observable of the transverse momentum, particularly 
the broadening of the transverse momentum distribution of Drell-Yan pairs produced 
on nuclei, as compared to a nucleon target. Empirically, the increase of the mean 
squared transverse momentum <p|.> is proportional to L^, the mean length of the 
quark path in a nucleus before creation of the Drell-Yan pair. 

5 <pl>j)Y= <pI>dy - <pI>dy= t^PAL (1) 

with K, being a dimensionless and A-independent constant; pA is the nuclear density. 
Relation (|l]) resembles broadening originating from a "random walk" -process. It is 
the purpose of our letter to derive eq. (|l|) and to calculate the constant k, starting 
from basic available information about strong interaction dynamics. 

Confinement is a fundamental aspect of strong interactions. The colour of a 
quark is always compensated by the colours of other accompanying partons. While 
this fact is unimportant at high momentum transfers, it is crucial for soft processes. 
Here the colour screening cuts off all low-momentum gluons whose wavelength is 
longer than the screening-length scale. Mathematically, this is equivalent to a cut- 
off of undesirable infrared divergences. 

We aim at calculating the transverse momentum distribution of a high-energy 
quark, which has traversed nuclear matter. We call this quark, originating from a 
parent high-energy hadron, the tagged quark in distinction to the screening partons 
which we call for simplicity "antiquarks" . At sufficiently high energies, the quark- 
antiquark separation is "frozen" by the Lorentz transformation, while the two par- 
tons pass the nucleus. These physical ideas can be technically realized as follows. 
The individual amplitudes of quark-nucleon scattering convolute to a hadron-nucleus 
multiple-interaction amplitude, with the tagged quark having in the final state an 
additional transverse momentum g, relative to the initial hadron momentum. The 
corresponding diagram is shown in fig. 1. Let us start with the n-fold scattering 
amplitude. Its square summed over all final states of nucleons participating in the 
interaction and integrated over all the transferred momenta, fixing the tagged-quark 
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final momentum g, has the form: 
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While n gluons are exchanged between the target nucleons and the incoming qq 
system, only ra (< n) gluons couple to the tagged quark (c.f. fig. 1). The wave 
function of the incoming hadron depends on the transverse momentum k of 

the tagged quark, and the relative longitudinal momentum < a < 1. The function 
/(fc) in eq. includes the two-gluon-nucleon coupling and gluon propagators 

m-\^0^^^-F.^^)\, (3) 

where F]^{k'^) is the two-quark formfactor of a nucleon: Fj^ikP') =< exp(ifc(fi — 
r2))\N>. The effective gluon mass irig is introduced to take into account confinement 
and is inessential at large k. We use the one-loop approximation for the QCD 
coupling constant as{k) at large fc, and freeze it at small k. The details can be 
found in 

After transformating the wave functions 0(A;|a) in eq. @) to the mixed p, a- 
representation, where p is a quark-antiquark separation in the impact parameter 
plane, eq. (H) changes into 
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da I d^p d^p ' e'^^P-p'^(t){p\a) 0*(p » [a{p, p')]" • (4) 



VAy„ (2vr) 
Here 

P') = / (0 /(^) (l - ^'''') (l - ^'') ■ (5) 
Note that o"(p, p ') can be represented as 

111 
o^(P,p') = 2^(^ + 2^^P'^ " 2^(P- P')^ (6) 

where 

is the inelastic interaction cross section of gg-pair with separation p, with a nucleon 

i- 

We sum eq. (^) over n with weight factors T"'{b)/n\, where T{b) is the nuclear 
thickness function at the impact parameter b, T{b) = pA{b, z) dz. In the Drell- 
Yan process one has to include the term with n = 0, which corresponds to no 
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initial state interaction of the incoming hadron in the nucleus before lepton-pair 
production. One also has to take into account the absorptive corrections, which in b- 
representation have a simple form, exp{ — ^[a{p) + a{p')] T{b)} being the probability 
amplitude to have no inelastic interactions in the initial and final state. 
Putting the factors together and using eq. (H) we get 
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The integration over q at fixed impact parameter gives 1, since (t(0) = 0, eq. (^. 
This result could have been anticipated since the quark does not disappear in the 
integration of all momenta. 

Expression eq. is the main technical result of this paper. It is of more general 
nature than the double-gluon approximation used as an input. One can use any 
QCD-inspired model for cr(p). Nevertheless we will use the results of the two-gluon 
approximation with running coupling as{k), of ref. [|l|. For small p^, <^ 1 fm^, 
where perturbative QCD is valid, one gets 

a{p) = Cp\ (9) 

with a dimensionless constant C = 3.2 The limit of small p is given in Drell-Yan 
production of heavy lepton pair to be discussed below, since the virtuality of the 
tagged quark is of the order of the squared mass of the pair. 

If in addition we use a Gaussian approximation for the hadron wave function in 



transverse direction 0(r) = y /i^/vr exp(— ^/i^r^), we obtain the cross section eq. (||) 
in the explicit form 

d^q J \[pl + 2CT{b)]'' • 

With the help of eq. (|10D, one can derive an expression for the cross section of the 
Drell-Yan process on a nucleus. We combine the transverse momentum distribution 
of the tagged projectile quark after it has traversed the nuclear thickness T(6, z) = 
Jl^dz' pA{b, z') with the transverse momentum distribution of an antiquark in a 
target nucleon to obtain 

^ = fd% HdzpAib^z) ^ , e-^V[<Pl>'h-y+2cnb,.)] _ 

d^q J J-oo ' n[<pl>%^ +2CT{b,z)] 

(11) 

Here, cr^y is the total cross section of the Drell-Yan process on a nucleon and 
<PT>Dy mean squared transverse momentum of the lepton pair produced in a 
hN collision. 



4 



One deduces from eq. ( [Tl| ) that the increase of the mean squared transverse 
momentum of the lepton pair produced on a nucleus compared to a nucleon as 
target, is 

5 <pl>),i.=<pl>),\ - <p|>^^= C <T>, (12) 

where < T >= ^ Jd'^b T'^{b), is the mean nuclear thickness. The result eq. ([T2|) 
confirms the empirical observation mentioned earlier eq. that the broadening of 
<Pt>dy of ^ lepton-pair produced on a nucleus is proportional to the mean length 
of path of the projectile quark in the nucleus, La =<T> /2pA- Moreover, 

K = 2C. (13) 

Eq. (|T3| relates the constant of proportionality k in eq. ([^) for the mean squared 
momentum transfer 6 <p'^> to the proportionality constant C, which governs the 
inelastic cross section a{p) = C eq. (|]) for the interaction of a qq pair transverse 
separation p with a nucleon. It is remarkable that according to eq. ([13D the constant 
K is independent of properties of the colliding hadrons and of the kinematics of the 
Drell-Yan process. 

The experiment NAlO |^ has measured a value 

5 <p|>^^= 0.15 ± .OSstat ± .03syst (GeV/c)^ (14) 

for incoming pions at 140 and 286 GeV, while the experiment E772 with 800 GeV 
protons IQ reports a value 

8 <pI>p^= 0.113 ± .016 {GeV/cf. (15) 

The two values coincide within the error bars which implies that no dependence on 
the type of incident hadron and its energy is visible, as expected from our expression 
eq. ([T3|). The expression eq. (0) predicts a value 

6<p^T >%= 0.17(GeV/c)' (16) 

for C = 3.2 from |I[ in fair agreement with the data. 

We calculate the ratio of differential cross sections of the Drell-Yan process on 
nucleus to nucleon targets, R{A/N), normalized with atomic weight A. We use 
expression ([TT|) with <Pt>dy from the systematics of of D-Y data. 

<p|>^P, = 0.59 + 0.0029 s 

< pj, >PPy = 0.39 + 0.028 yi. (17) 

Here, s is the square of the total cm. energy in GeV^. Assuming a Gaussian form 
for the pr-distribution, we use <p'^>= ^ < pr 
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The results of calculations of R{A/N) are compared with experimental data 
[0, on fig. 2. One can see that our parameter-free calculations provide a good 
description of experimental data. The flattening of the p^-dependence of R{A/N) 
at high energies observed experimentally is a result of the increase of <p^>^^ as a 
function of energy. 

Conclusions 

• The phenomenon of broadening of the transverse momentum distribution of 
a quark propagating through nuclear matter crucially depends on the effect 
of colour screening in soft interactions. While the cross section cr(p) for the 
individual quark-nucleon scattering decreases with decreasing screening ra- 
dius p, the mean transverse momentum picked up in each collision is inversely 
proportional to p. The product of <Pj'> ■'^{p) in each collision is therefore in- 
dependent of p and equals the constant 2C. The smaller the screening radius, 
the rarer the quark rescatters, but the larger is mean transverse momentum 
in each rescattering. 

• We calculate the cross section a{p) in the one-gluon exchange approxima- 
tion for the hadron inelastic scattering amplitude. This approach is justified 
because of the high virtuality of the quarks participating in the Drell-Yan 
reaction. 

• Analogous consideration may be applied to the hadroproduction of charmonia 
on nuclei which predominantly proceed via gluon fusion. We expect an ap- 
proximate doubling of the broadening effect due to the additional factor 9/4 
in the cross section of the gluon-nucleon interaction. Experimental data on 
J/\l' production on nuclei at 200 GeV confirm the expected increase and give 
values for kj/^i/kdy = 1-9 ± 0.4 and 1.6 ± 0.25 [0]. However, for the T 
production at 800 GeV 0] the experimental ratio kj/^/kdy = 5.9 ± 1.3 is 
much larger than what we would expect and needs another explanation. 

• The process of inclusive production of hadrons with high transverse momenta 
off nuclei ("Cronin effect") has much to do with the colour screening effect. 
Indeed the high transverse momentum is built up of many soft scatterings in 
a way similar to the case treated in this paper for the Drell-Yan process. 

• The broadening of the transverse momentum of a quark propagating through 
nuclear matter, might be relevant also to the process of inclusive production 
of leading hadrons in deep inelastic scattering on nuclei. However, it occures 
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only if the hadron is produced outside the nucleus, i.e. the length of the 
formation zone exceeds the nuclear radius. Otherwise, the broadening becomes 
much weaker than in the Drell-Yan process, as the colourless hadron should 
reinteract in the nucleus with a small elastic cross section. 
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Figure captions 

Fig. 1: The inelastic interaction of a quark q and a antiquark q with n nucleons of 
target nucleus. 

Fig. 2: The ratios R{A/B) of the Drell-Yan cross sections da^ / d'^px and da^ j d'^px 
as a function of the transverse momentum pt of the D-Y pair. The data from 
experiments with various projectiles, targets and energies are compared with the 
calculation based on eq. (p4]). 
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